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Abstract. We show that the short-scale elastic distortions that are excited in the vicinity of the joint be- 
tween diflterent lipidic membrane domains (at a scale of ~ 10 nm) may produce a "crease" from the point of 
view of the standard elastic description of membranes, i.e., an effective discontinuity in the membrane slope 
at the level of Helfrich's theory. This "discontinuity" may be accounted for by introducing a line tension 
with an effective angular dependence. We show that domains bearing strong spontaneous curvatures, such 
as biological rafts, should exhibit creases with a finite contact-angle, almost prescribed, corresponding to 
a steep extremum of the line energy. Finite contact-angles might also occur in symmetric membranes from 
the recruitment of impurities at the boundary. 

PACS. 87.16. Dg Membranes, bilayers, and vesicles - 87.10.-|-e Biological and medical physics: general 
theory and mathematical aspects - 68.03.Cd Surface tension and related phenomena 



I 1 Introduction 

I When phospholipid or like surfactants are dissolved in an 
' aqueous solution, they condensate into fluid bilayer mem- 
, branes. Lipid bilayers are formed by two contacting mono- 
■ layers of opposite orientations, in which the hydrophilic 
I heads of the surfactants are directed towards the aqueous 
. solution while the hydrophobic tails, confined within the 
' sheet, are shielded from contact with water [1]. Lipid bi- 
, layer membranes are ubiquitous in biological cells ,2 . In 
' actual biological systems, membranes are formed of multi- 
' pie lipid or surfactant components. The latter may later- 
, ally separate into coexisting liquid phases, or domains [3', 
|2n5j , with distinct compositions and distinct microscopic 
] features, such as the membrane thickness. Microdomains 
called "rafts" are receiving increasing attention, since they 
are believed to concentrate proteins that must interact 
with one another to carry out important cellular func- 
tions [6l[7]. 

A number of theoretical works have studied the re- 
lationship between the curvature elasticity and line ten- 
sion of domains and the resulting formation of particular 
patterns or membrane shapes [8MfT0l[TTl[T2lfT3lfT4] . In all 
these works, the slope of the membrane is assumed to be 
continuous at the domain boundary. In this paper, we ar- 
gue that although this assumption is realistic in a large 
number of situations, it is not correct for domains with 
strongly asymmetric monolayers, as it is often the case in 
the biological realm [T5] . 



Here, we need to precise what we mean by a mem- 
brane "slope discontinuity" : we mean a discontinuity as 
far as the standard Hclfrich description of membranes is 
concerned. This "discontinuity" will obviously be resolved 
at a more microscopic scale, i.e., at the molecular scale or 
at the scale of more involved elastic theories taking into 
account short-scale structural degrees of freedom. We do 
not refer here to apparent large-scale discontinuities as, 
e.g., the "neck" connecting a budding domain in a vesi- 
cle f9l[T0]: these are actually smooth and continuous at a 
sufficiently short scale within Helfrich's theory. The possi- 
bility of actual slope discontinuities in Helfrich's theory, or 
"creases" , was suggested on symmetry grounds for surfac- 
tant monolayers in Ref. [16| . but not investigated further. 
Here, we actually study the possibility of such creases in 
the case of membranes. 

At first sight, the condition of a continuous slope across 
the domain boundary seems mathematically necessary be- 
cause of the presence of a curvature energy. Indeed, in the 
standard Helfrich model [17], which describes the mem- 
brane as a fluid structureless surface with a curvature 
elasticity, the energy density is proportional to the square 
of the membrane curvature. (If the membrane is asym- 
metric, there also a term linear in the curvature, with a 
coefficient proportional to the asymmetry of the bilayer.) 
Within this simple model, which is essentially valid at 
large length-scales, if one assumes a discontinuous mem- 
brane slope at the boundary between two domains, one ob- 
tains a localized infinite energy (the local energy density 
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diverges there as the square of a delta- function) . However, 
this reasoning assumes that the Helfrich model may be 
continuously applied through the boundary between the 
domains, which is not obvious. In fact, as in every elas- 
tic model, there is a cutoff below which Helfrich's model 
is no longer applicable. The corresponding length-scale is 
comparable with a few times the membrane thickness, i.e., 
~ 10 nm [18. . At the joint between membrane domains, 
there are short-scale elastic degrees of freedom that are 
excited: matching of the two membrane thicknesses |19| 
[5D], tilt of the lipidic tails 21j, etc. These degrees of free- 
dom relax typically on a length scale comparable with the 
cutoff of the Helfrich model. In this paper, we show that 
these structural degrees of freedom may result in an ap- 
parent contact- angle in the large scale Helfrich description 
(an effective jump in the membrane slope at the joint). We 
also demonstrate that the energy of these extra degrees of 
freedom may be "coarse-grained" into an anisotropic line 
energy 7(0) for the Helfrich model, where 9 is the effective 
slope discontinuity at the joint. 

2 Model 

Let us consider two contacting membrane domains Vi, 
i G {1,2}. We shall describe them more "microscopically" 
than in Helfrich's theory. Since we are interested in the 
local junction between these domains, we assume trans- 
lational invariance along the y-direction. We take Vi = 
[— L,0] and I?2 = [0, L], the joint being at a: = 0. Our 
model free-energy, per unit length, is the following: 




The first contribution /jj is similar to the bending energy 
density of the Helfrich model [I7j : 

fl^^ ^l^.{d^h,f -K,c,d'h,. (2) 

Here hi{x) represents the height profile of the membrane 
midplane in domain i. Dealing only with small deforma- 
tions, we shall neglect everywhere terms of higher-order 
than hf. With d = d/dx, the quantity d^hi represents the 
membrane curvature, and Ki and Ci are the bare bending 
rigidities and bare spontaneous curvatures, respectively. 
The latter term, which is linear in hi, is allowed by sym- 
metry only if the membrane is dissymmetric, i.e., if the 
lipid or surfactant compositions in the two monolayers 
are different. For the sake of simplicity, we neglect the 
surface tensions ai , which would give rise to a term of the 
form ^ai{dhi)^: this is reasonable as long as we investi- 
gate length scales smaller that (Kij Oi)^!"^ . Note that the 
Gaussian rigidities Ri [T^ can be discarded, since in one 
dimension the Gaussian curvature vanishes. 

The second term, , arises from the structural en- 
ergy density associated with the inner deformations of the 
membranes. Although the length-scales under consider- 
ation are comparable with a few nanometer, we keep a 



continuous description. Among the many possible struc- 
tural variables, two traditional ones are the thickness of 
the membrane and the tilt of the lipids relative to the 
membrane's midplane (see, e.g., pT|[22]). Note that our 
aim is not to discuss the most general model, but to show, 
with the simplest symmetry-allowed terms, the possible 
existence of the "contact-angle" effect discussed in the in- 
troduction. Hence, for the sake of simplicity, we consider 
only the membrane thickness variable: 

= ^B, {u, - a,f + \h [du^f . (3) 

Here, Ui(x) describes the thickness of the membrane in 
domain i. These two terms represent the free-energy cost 
associated with a variation of Ui with respect to its equi- 
librium value a^, at lowest-order in the gradient expan- 
sion [22ll23il24j . We expect the typical length-scale {ki/BiY^^ 
over which Ui relaxes, i.e., the width of the joint, to be 
comparable with the membrane thickness. 

There could be no possible interplay between the struc- 
tural degrees of freedom and the large-scale shape of the 
membrane without an energy term coupling these quan- 
tities. We consider the following, lowest-order interaction 
term: 

fl^^-A,{u,-a,)d^h,. (4) 

It represents the coupling between the excess thickness 
and the spontaneous curvature. Like the spontaneous cur- 

(i) 

vature term in , this term is allowed by symmetry only 
if the two monolayers forming the membrane are dissym- 
metric. 

Finally, 70 is a line tension term that takes into ac- 
count the microscopic interaction between the two do- 
mains (arising, e.g., from van der Waals forces). 

For the sake of simplicity, we assume Bi = B2 = B, 
fci = fc2 = fc and A\ = Ai = A, but we keep ki 7^ K2, 
ci 7^ C2, and ai ^ 02- We shall discuss later on the or- 
ders of magnitude of the different terms. Note that the 
two-dimensional version of this model would be simply 
obtained by replacing dx by d^x and 9 by V. 



2.1 Renormalized bending rigidity and spontaneous 
curvature 

A straightforward effect of the coupling term A{ui—ai)d'^hi 
is to renormalize, away from the joint, the bending rigidi- 
ties and spontaneous curvatures |25]. Indeed, if we ne- 
glect dui (by assuming that the thickness is uniform far 
from the joint), the total free-energy density reduces to 
\K,i{d'^h^Y - Kidd^hi + \B{ui - - K^i - ai)d'^hi. 
The latter yields, after minimization with respect to u^, 
Ui = ai + {A/ B)d^hi, and the energy density becomes 
^K'^{d^h,f - K'^c'^d^h^, with 
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These are the effective bending rigidities and spontaneous 
curvatures of the membrane : the ones that would be 
measured in a macroscopic experiment and the ones that 
should appear in the Helfrich curvature energy. 



2.2 Equilibrium equations 

The Euler-Lagrange equilibrium equations associated with 
our elastic free-energy F are 



Ki d^hi{x) 
k d'^Ui{x) 



Ad'^Ui{x) , 
B [ui{x) - a.-] 



Ad'^hiix). 



(7) 

(8) 



These equations should be solved with the correct number 
of boundary equations. First, there are three continuity 
equations at the joint: 



{h2 
{h2 



\x=Q 
\x=Q 



{hi 
{hi 



\x=0 
U=0 



dh2\x=o = dhi\^=Q. 



(9) 
(10) 
(11) 



The first two conditions represent the requirement that 
the polar heads of the two membranes must match at the 
joint. They imply the continuity of both hi and Ui. The 
third one, the slope continuity equation, must be imposed 
in order to avoid an infinite curvature energy, as explained 
in the introduction. Note that there is therefore no contact 
angle at the "microscopic" level in this model. 

When calculating the variation SF associated with the 
variations Sui{x) and Shi(x), there are boundary terms in 
x = in factor of 6u\x=q, Sh\x=Q and Sdh\,j;=Q (remember 
that u, h, and dh are continuous at the joint). These are 
boundary forces and torques. The associated equilibrium 
conditions are 



fl\x=0 + f2\x=0 — , 
fl\x=0 + f2\x=0 = , 
i^l\x=0 + ^^2\x=0 = , 



(12) 

(13) 
(14) 



corresponding to the vanishing of the total force relative 
to h, the total force relative to u, and the total torque 
relative to h, respectively, with 



fi = -e {Ki d^h, 
fl^ ek dui , 
fii = e [ni {d'^hi - 



A dui 



i) ~ A {ui - a.i)\ 



(15) 
(16) 
(17) 



Here, i £ {1,2} indicates the domain where the quantity 
is calculated, and e = 1 if the quantity is calculated at the 
right-hand side boundary of its domain, while e = — 1 if it 
is calculated at the left-hand side. 

In addition, we assume that the domains "Di actually 
correspond to small regions in contact with larger mem- 
brane parts, in equilibrium under some unspecified bound- 
ary conditions and global constraints. These outer parts 
transmit boundary forces f^^^\x=±L, U=±l and tor- 



ques l?o 



^±L, which yield the equilibrium equations: 



f2\x=L + /° 



0, 



(18) 



fl\x=-L + P 



f2\x=L + /'™'U=L 



0: 

f[\x = -L + r''X = -L^Q., 
^2\x=L + ^cxt\x=L = : 



(19) 
(20) 
(21) 
(22) 
(23) 



These conditions can be formally derived by adding to F 
the contribution : 

Fext = n = -L h{-L) + r%^L h{L) 



-_^Ldh{-L) + 



cxt \x—L 



dh{L), (24) 



and by equilibrating the boundary terms. Notice that Eqs. 

correctly provide 12 boundary conditions for two 
integration domains subject to differential equations glob- 
ally of the sixth order. 



2.3 Equilibrium solutions 

The general solutions of the set of linear differential bulk 
Eqs. Q-® are: 



hi = A,i + Ai2 X + -Aa x^ + -An x^ 



A. 



i5 2 



e"' +A, 



A 



(25) 



u, ^a, + ^Aa + '^A,iX + Ai5 ^ + A,e e^^" ^, (26) 
where Aia are integration constants, and 



(27) 




ii ^ \ — \ -r 



Since qi w ^ B/k, we expect it to be of the order of the 
inverse of the membrane thickness (unless when — > 0). 

We see that the solutions are the sum of polynomi- 
als, which correspond to the solutions one would get in 
a pure Helfrich model plus relaxing exponentials, which 
correspond to the effect of the structural degrees of free- 
dom. Since these exponentials appear also in hi, they may 
affect the extrapolated angle at which the polynomials so- 
lutions meet the joint. This is the source of the effect we 
are interested in. 



3 Detailed structure of the joint 

To study the structure of the joint and its generic response 
to the constraints transmitted by the curvature elastic- 
ity of the membrane, it is enough to consider the situa- 
tion in which two opposite far-away torques are applied: 
— = Q. Remember that those torques 

are acting at the boundaries of the domains Di and that 
they represent the effect of the pieces of membranes that 
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are outside the region we are studying. We set to zero 
the external boundary forces: f°^*\±L — f'°^^\±L = (the 
latter would merely transmit extra torques). We also as- 
sume that the membrane domains are much larger than 
the "width" of the joint, i.e., L ^ q^^. 

3.1 Solution of the problem 



With the above conditions, it follows from Eqs. 
that Aj4 = 0, AiQ = A25 — and that 



A,: 



(28) 



Obviously, we may choose A21 — A22 = by a simple 
translation and rotation of the reference frame. Requiring 
that the remaining boundary conditions be satisfied yields 



with 



All = 

A12 = 
Ai, - 

^26 - 



A 

a2 + - 



A 



K2 ql 



Kiqf 



«iK2 qfq2 {qi +q2) 
Afi B{n2-Ki) 



K1K2 kqi q2 {qi 
92 



92) 



qi - 
qi 



92 



qi + 92 



(29) 
(30) 

(31) 
(32) 

(33) 



We therefore find that the membrane shape hi relaxes 
exponentially away from the joint with a characteristic 
length q~^ to a parabola of curvature Ais, equal to the 
effective spontaneous curvature plus a deviation pro- 
portional to the applied external torque (Eq. (pS)) ). The 
membrane thickness, Ui, relaxes exponentially to the value 
ai + {A/B)Ai3, which is shifted with respect to Oj by an 
amount proportional to the curvature Ai^ . These features 
can be seen in Fig. [1] Away from the joint, hi{x) tends to 
Hi{x) = Ail + Ai2 X + |^i3 2^^, as in the Helfrich model 
(for small membrane inclinations). 

Although the joint essentially appears as a thickness 
matching, by extrapolating the large-scale solutions Hi{x) 
up to x = 0, one finds a non-zero contact-angle (Fig. [T]), 
e = 9i/iU=o - dH2\x=Q, given by 



= Ai2^0. 



(34) 



As expected, 6 is proportional to A: the coupling between 
the membrane shape and the structural degrees of freedom 
is necessary to obtain an effective contact-angle. Notice 
that there is also a mismatch in the extrapolated mem- 
brane height, Hi\x=o — H2\x=o 7^ 0, but we may neglect it 
because it is smaller than the membrane thickness. 



3.2 Energy of the equilibrium configuration 

By integrating F by parts and making use of the bulk 
equations, the total energy at equilibrium, F = Fq, may 




X [nm] 

T2l) - (|23p Fig. 1. Structure of the membrane in the vicinity of the joint. 

The upper and lower bold lines correspond to hi{x) + ^Ui{x) 
and hi{x) — ^Ui{x), respectively, where i = 1 for a; < and 



i = 2 for X > 0. The parameters are K2 = 2x 10 
ci = (50nm)~\ C2 = 0.03 ci, B = 4x 10^^ Jm~*, ai = 7nm, 
a2 = 5nm, k = 10"Mm~^ A=1.2x 10"^ Jm~^ and i7 = 0. 
The asymptotic large-scale solutions Hi{x) are displayed as 
thin lines: their extrapolated contact-angle is 9 = 9o — 17°. 
For a better view, the membrane has been rotated hy 0/2 coun- 
terclockwisely (and Ui has been corrected in order to yield the 
correct thickness when the slope is large). 



be expressed as boundary terms only: 
-^0= l[f[iui-ai) + fihi + if2i-KiCi)dhi] 



x=0- 



X! 9 [■^2 ("2 - 0,2) + /2 h2 + {^2 - K2 C2) (9/l2] 



(35) 

Using the boundary conditions ((9|- (|23|) . and substituting 
the solution while taking the limit Lqi 3> 1, the total 
energy may be cast into the following form: 



70 + ll& + ]p2 I 



where 



71 



72 = 



^{ai 



0-2) 




(36) 

(37) 
(38) 



and 7o is the constant initially present in F. One recog- 
nizes in Eq. ([36]) the effective (large-scale) Helfrich energy 
calculated from the asymptotic curvatures A^i (which are 
constant here) and the renormalized bending rigidities and 
spontaneous curvatures k[ and , plus an effective line en- 
ergy 7(6*) = 716* 1726*2. 

Note that the large-scale curvatures A^i^s and the ef- 
fective contact-angle 6 are controlled by the applied torque 
n (see Eqs. ([Mll-dMl)). In fact, from Eqs. ^ and ®, 
one can show with a little algebra that 



n = --f2ie - Oa) , 

where 9q is the extremum angle of 7(0), given by 



(39) 
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ai-a2 + — (ci - C2) 



71 

72 



(40) 



Eq. simply shows that the external torque f2 (which 
is transmitted up to the joint by the bending rigidity of 
the membrane) is equilibrated by the effective torque gen- 
erated by the boundary line. 

Note that we obtain 72 < 0. Because of this sign, the 
function "f{0) possesses a maximum in 6 = 0o and not a 
minimum. This is not a problem, actually, because 7(0) is 
only an excess energy, which does not exist by itself in the 
absence of a surrounding membrane. Hence, there is no 
stability criterion for 7(6*) alone [15]. The only condition 
physically required is that the total energy be minimum 
for i? = 0, or, equivalently (since 9 is an affine function of 
n) for 9 = 9q. This is indeed the case for L > q~^, since 
the first term of Eq. (j36p , which is obviously minimum for 
— Ci, i.e., for Q ~ Q, i.e., for 9 = is proportional to 
L, while 7(0) is not. 



4 Effective large-scale Helfrich model with a 
boundary crease 

We have therefore shown that our "microscopic" model 
taking into account the membrane thickness structural de- 
gree of freedom is equivalent, at the coarse-grained level, 
to the following (pure) Helfrich model supplemented by an 
anisotropic line energy for the effective joint contact-angle 
9 = dhi\ 

x=o ~ dh2\x=o (that must be considered as free to 
adjust to equilibrium): 



dx 



+ 7o + (^ - ^0) 



(41) 



where 7g = 70 — ^71/72 > 0. Here, hi{x) does not describe 
any longer the height profile of the membrane midplane, 
it simply describes the membrane as a whole, since at the 
level of Helfrich's description the whole bilayer is treated 
as a pure mathematical surface. 

Let us detail the equivalence between the two models 
at the coarse-grained level. The equilibrium equations and 
boundary conditions at the joint, for the effective Helfrich 
model with energy are 

d^h,{x)=0, (42) 

h2\x=o = hi\x=o , (43) 

dh2\x=o = dhilx^o + 8 , (44) 

/iU=o + /2U=o = 0, (45) 

^iU=o = -7'(^)--^2U=o, (46) 

with, now, fi — —eKid^hi and f2i — eKi(d^hi — Ci), e be- 
ing defined as previously. We see indeed that since dui 
tends to zero far from the joint, Eq. ([7]) is equivalent at 
the coarse-grained level to Eq. (|^^ . and that the effective 
torque equilibrium equation (|39p is equivalent to Eq. (|46p . 



The continuity Eq. (|44|1 is postulated by the very def- 
inition of 9 (which becomes a free parameter). Finally, 
Eqs. (|43l) and ([45]) are only approximately verified, but 
with a very good precision since we have seen that the 
height mismatch in the extrapolated large-scale solutions 
is smaller than the membrane thickness. 

Note that Eq. (jlS)) shows that 9^ is the equilibrium (ef- 
fective) contact-angle of the joint, attained when no torque 
is transmitted by the membrane, i.e., when the membrane 
is curved on both sides of the joints at its spontaneous cur- 
vature. This corresponds to the situation shown in Fig. [T] 

In the limit yl ^ (in which Ui becomes decoupled 
from the membrane shape hi), we should recover the stan- 
dard slope continuity condition. This is indeed the case, 
since Eqs. ([38]) and ([40]) yield then ^ and 72 ^ 00. 
We also obtain 7o ^ 7o + jVBk (oi — 02) , which is sim- 
ply the sum of the bare line tension and the elastic energy 
associated with the thickness mismatch in the absence of 
any coupling. 

Finally, note that we have verified the equivalence of 
the large-scale behavior of models F and F'^^ for arbi- 
trary external forcings (including external forces in addi- 
tion to external torques) by numerically solving the two 
boundary-value problem and comparing the solutions far 
from the joint. 



5 Orders of magnitude 

Let us discuss the orders of magnitude of the various pa- 
rameters for a joint between two domains with a large 
thickness mismatch, e.g., as it is the case in the junction 
between biological membranes and rafts [T9ll20j . For bio- 
logical membranes K2 = 2 x 10~^^J and 02 = 5nm are 
standard values [T5 . Since rafts are known to be thicker 
and more ordered than normal membranes, we take ai — 
6nm and ki = 5 K2- From the area-stretching coefficient 
kg ~ 0.1Jm~2 [X] and the volume conservation relation 
B ~ fcs/a|, we obtain B ~ 4 x 10"^^ Jm~'*. To estimate k, 
we assume that the relaxation length of Ui is comparable 
with the membrane thickness, which yields k kg. These 
values are used in Fig. [T] 

If the two monolayers are identical, we must set = 
and A = Q hy virtue of the resulting up-down symme- 
try, and hence there is no contact-angle effect within the 
present model. For asymmetric membranes, the values of 
the spontaneous curvature radii may range between 
several hundred of micrometers to a few tenth of nanome- 
ters, depending on the difference in lipid composition or 
inclusion concentrations between the two monolayers. 

We must now estimate the parameter A. By compar- 
ing the terms Ki Ci&^hi of /h and A{ui — ai)d^hi of fmt, we 
see that A measures the dependence of the spontaneous 
curvature in the membrane thickness. We expect this de- 
pendence to be significant, since it is well-known that the 
spontaneous curvature of monolayers is (qualitatively) re- 
lated with the effective conical shape of its constituents [1] . 
In the absence of any other information, let us assume 
that a 10 % change in the membrane thickness a (which is 
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Table 1. Values of tlie asymmetry-dependent constant A and 
the corresponding joint parameters as a function of tlie mem- 
brane spontaneous curvature radius c~^. 



c ^ [nm] 


A [pN/nm] 


1^0 1 [deg] 


70 [PN] 


R* [nm] 


50 


20 


23 


70+5.002 


96 


300 


3.3 


2 


70+5.001 


0.7 


10* 


0.01 


6 X 10~^ 


70+5 


6 X 10"^ 











70+5 






a large variation) may change the spontaneous curvature c 
by an amount comparable with the spontaneous curvature 
itself. This yields the estimate yl x 0.1 a ~ kc, which, with 
K ~ 5 X 10~^^ J and a ~ 5 nm, gives A ~ 10^^ J m^^ x c. 
As shown in Table [TJ it follows that significant contact- 
angles should only be expected for strong spontaneous 
curvatures. We see also that the contribution to the ef- 
fective line tension depends very weakly on A. 

5.1 Rigidity of the crease 

To determine the angular rigidity of the joint, we can es- 
timate the curvature radius R* which must be applied to 
the joint in order to obtain an angular variation of 0.1 rad 
(~ 5.7 deg). From Eq. we obtain R* ~ K/{-f2 x .1) 

(see Table [J) . We see that even for large values of A, R* is 
quite microscopic; hence, the joint appears as extremely 
rigid. In consequence, we may almost view Oq as a pre- 
scribed contact- angle. 



6 Discussion and conclusions 

The present study shows that if the structural distortions 
arising within the joint between membrane domains are 
coupled with the membrane shape, then one should gener- 
ically expect, in the large-scale elastic description, an ef- 
fective slope discontinuity 9. The boundary line energy 
should then be considered, in general, a function ^{9). 

It should be noted, however, that if an extremely sharp 
extremum occurs in 7(0) for = 0, then one effectively re- 
covers the usual slope continuity condition and one should 
not bother about the angular dependence of j{9). The sit- 
uation is therefore subtle and it is essential to consider the 
conditions for which either 9q or the extremum of j{9) 
is shallow (compared to the torques produced by macro- 
scopic elastic distortions). 

It is not likely that the extremum of 7(6*) be effec- 
tively shallow. Indeed, finite apparent angular discontinu- 
ities must be resolved within the joint by strong contin- 
uous distortions, and the latter will always cost energies 
very large compared to those needed to curve the mem- 
brane on a macroscopic scale. It is thus more promising 
to look for ^0 7^ 0, but this is possible only for membranes 
with broken up-down symmetry, hence membranes bear- 
ing a spontaneous curvature. This is one of the reasons for 
the choice of the specific structural model studied here. 



In fact, in this paper, we have examined the conse- 
quences of the most elementary boundary distortion: the 
membrane thickness matching (occurring, e.g., at a raft 
boundary). For the sake of simplicity, we have consid- 
ered homogeneous membranes and we have assumed a 
linear coupling between the thickness and the shape of 
the membrane (the latter requiring the broken up-down 
symmetry). We have demonstrated, within this model, 
the existence of an effective anisotropic line tension 7(6*) 
and we have shown that it has such a sharp extremum 
^0 that the latter can be considered almost as a pre- 
scribed contact- angle. We have obtained 6*0 — for mem- 
branes with macroscopic curvature radii, but finite, for 
membranes with microscopic curvature radii. This is in- 
deed interesting, because biological membranes, and espe- 
cially rafts are known to be very asymmetric, with sphin- 
gomyelin and glycosphingolipids enriched in the exoplas- 
mic leaflet and glycerolipids in the cytoplasmic leaflet [6] . 
Curvature radii of a few tenth of nanometers, as in TableO 
are thus plausible, and biological rafts could be expected 
to have contact-angles as large as 10 — 20°. 

One should bear in mind that we have used a contin- 
uous model to investigate the effect of nanometric scale 
distortions. This is actually the limit of validity of con- 
tinuous membrane models, and one should not expect 
better than semi-quantitative results. Besides, there are 
many other non-linear terms that might be taken into 
account, and also many other structural degrees of free- 
dom. Real contact-angle might therefore be significantly 
stronger, or smaller, than what we have estimated here; 
experiments, or numerical simulations (e.g., with coarse- 
grained molecules) should be the best way to give a quan- 
titative answer. Note also that we have neglected, for the 
sake of simplicity, the higher-order terms proportional to 
d^u and (d^u)^ which have been investigated in [21]: they 
certainly would renormalize our results but cannot, for 
symmetry reasons, be the direct source of an effective 
contact-angle. 

An interesting effect that we have neglected here is the 
possibility that the boundary will recruit dilute impurities 
and accumulate them in the joint, in order to relax some 
of the microscopic distortions. Since such impurities will 
in general have different concentrations in the two mono- 
layers, they will locally enhance, within the joint, the local 
spontaneous curvatures and thus the local coupling con- 
stant A, yielding again finite contact-angles (even for up- 
down symmetric host membranes) . These impurities could 
also soften the elasticity of the boundary distortions. 
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